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Abstract 

The shadowing of cosmic ray primaries by the the moon and sun was observed by the 
MINOS far detector at a depth of 2070 mwe using 83.54 million cosmic ray muons accu- 
mulated over 1857.91 live-days. The shadow of the moon was detected at the 5.6 a level 
and the shadow of the sun at the 3.8 a level using a log-likelihood search in celestial coor- 
dinates. The moon shadow was used to quantify the absolute astrophysical pointing of the 
detector to be 0.17±0.12°. Hints of Interplanetary Magnetic Field effects were observed in 
both the sun and moon shadow. 



1 Introduction & Motivation 



The Main Injector Neutrino Oscillation Search (MINOS) far detector 0]| is a magnetized scintil- 
lator and steel tracking calorimeter, located in the Soudan Mine, ( 47° 49' 13.3" N, 92° 14' 28.5" 
W) in northern Minnesota, USA, at a depth of 2070 meters water equivalent (mwe). While the 
primary function of the far detector is to detect neutrinos from Fermilab's NuMI beam [2J, the 
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great depth and wide acceptance of the detector combined with the fiat overburden of the Soudan 
site allow it to serve as an efficient cosmic-ray muon detector. The 5.4kton detector is composed 
of 486 8 m wide octagonal planes, each consisting of a 2.54 cm thick steel plate and a 1 .0 cm thick 
scintillator plane. Adjacent planes are separated by a 2.4 cm air gap. It is 30 m long and has a total 
aperture of 6.91xl0 6 cm 2 sr fl3) for the cosmic rays selected in this analysis. MINOS observes 
underground muons with a minimum surface energy of 0.7 TeV, and the sharply peaked energy 
spectrum has a mean value of about 1.0 TeV [4J. This mean muon energy corresponds to a mean 
cosmic ray primary energy of about 10 TeV. 

Optical telescopes use a standard catalog of stars to establish the resolution and pointing reliability 
of a new instrument. This is not possible for a cosmic ray detector as there are no known cosmic ray 
sources available for calibration [5|. The shadow caused by the absorption of cosmic rays by the 
moon is a well observed phenomenon in the otherwise isotropic cosmic ray sky [6|. This shadow 
provides a means of studying the resolution and alignment of the detector which are important in 
the search for cosmic point sources. The physical extent and shape of the shadow gives information 
about the resolution of the detector, while the location of the deficit center measures the absolute 
pointing of the detector. The moon has a 0.52° diameter as viewed from earth, and the cosmic 
ray deficit it causes has been measured by air shower arrays (CYGNUS [7], CASA 0, Tibet [|9]|, 
Milagro OH, GRAPES 03), HEGRA |23) as well as underground detectors (Soudan 2 [13], 
MACRO EMI, L3+C El, BUST ifTTll). 

The shadow of the moon is affected by multiple Coulomb scattering, the geomagnetic field and the 
Interplanetary Magnetic Field (IMF) [fT3l . Multiple Coulomb scattering occurs in the rock over- 
burden and causes a general spreading of the moon deficit disc. The geomagnetic field is nearly a 
dipole and causes an eastward deflection of positive primaries, which results in a positive horizon- 
tal shift in the observed shadow of up to Aa=0.15° Z j E p (TeV) (see Sec. I2.2I ). Older calculations 
using an impulse approximation give a larger value [1 1 81 1911 . The IMF is produced by the sun, which 
has an ambient dipole field that is 100 times greater than the geomagnetic field. The field is carried 
through the solar system by the solar wind, the stream of energetic charged particles that emanate 
from the atmosphere of the sun. Since the sun has a 27 day rotation period, the magnetic field has 
a spiral shape, called a Parker spiral [20] . The IMF causes a deflection of primaries that strongly 
depends on the solar wind. It changes in time and has a sectorized structure ED . This complex 
structure makes it hard to model. The direction of the IMF vector is either "towards" or "away" 
from the earth and the two field directions are separated by a thin, field-free, region known as the 
"neutral current sheet". The IMF causes a deflection that smears the moon's shadow, though this 
effect is small since a cosmic ray primary travels a relatively short distance from the moon to earth. 
The larger IMF effect comes from its distortion of the geomagnetic field. 

The sun also subtends a disk of diameter 0.52°, and in principle produces a similar shadow to the 
moon. Its shadow has been observed by CYGNUS Q, Tibet Il9l2"2l and MACRO (T5|. However, 
the much greater distance that cosmic primaries travel in the IMF can produce significant distor- 
tions in the shadow. The IMF varies according to the 1 1 year solar activity cycle and peaks at solar 
maximum, when the sun's magnetic field changes polarity [|23l . The most recent solar maximum 
occurred in February, 2001, and the following minimum occurred in December, 2008. The average 
magnitude of the IMF over the period of time that the MINOS data were collected was 6.5 nT. The 
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maximum value was 50 nT. Assuming a constant, uniform field and a lOTeV proton, the average 
magnetic deflection is about 0.05°. For a lOTeV proton traveling along the line from the Sun to 
the earth, the deflection is up to 0.75°. Because of the complex structure of the IMF already men- 
tioned, the assumption of a constant, uniform field is not correct. By measuring the behavior of an 
ensemble of particles, thus sampling many parts of the IMF, a more refined understanding of the 
average IMF properties can be obtained. 

1.1 Spherical Coordinate Systems 

Horizon coordinates use the detector's local horizon as the fundamental plane. The longitudinal 
angle is azimuth (Az), measured east from north and ranges from 0° to 360°. The angle out of the 
equatorial plane is zenith (Zen), where vertical up has a value of 0° and the horizon has a value 
of 90°. Zenith is measured from 0° to 180°. Celestial coordinates are projected on the sky and 
centered on the earth's equatorial plane. The longitudinal angle is right ascension (RA), measured 
eastward from the Vernal equinox and ranges from 0° to 360°. The angle out of the equatorial 
plane is declination (Dec), ranging between ±90°. Ecliptic coordinates are celestial coordinates 
that use the ecliptic (the path the sun follows over the course of a year) as its fundamental plane. 
The longitudinal angle is ecliptic longitude (A), and is the same as RA. The angle out of the plane 
of the ecliptic is ecliptic latitude (J3), ranging between ±90°. 



2 Data 

In order to perform the shadowing analysis, the muon data must be selected for reliable pointing, 
the smearing and systematic offsets of the shadow must be understood, and the backgrounds must 
be quantified. 

2.1 Event Selection 

This analysis encompassed events recorded over 1980 days, from August 1, 2003 to Decem- 
ber 31, 2008, giving a total of 1857.91 live-days when the detector was operational. The data 
set includes 83.54 million cosmic ray induced muon tracks. Cosmic ray muons were triggered by 
recording hits on four planes within a group of five planes. Several selection criteria were required 
to ensure that the detector was in a reliable state when the data were taken (Pre- Analysis selection 
criteria) and that only well reconstructed tracks were included in the sample (Analysis selection 
criteria). The Pre- Analysis selection criteria are described in Il24l . The following analysis selection 
criteriawere applied: 

(1) "Number of Planes > 9", a track that passes fewer planes may not give reliable localization 
information to the track fitter. 
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(2) "Track Length > 1.55 m", any event with a track shorter than 1.55 m may not be reliably 
reconstructed. 

(3) "\~o v tx — ~o end\ < 0.021m" If the uncertainty on the endpoint position, a en d, is signif- 
icantly different from the uncertainty on the beginning position, a vtx , then the muon has 
questionable reconstruction pointing. \a vtx — ~^ en d\ is a measure of the absolute value of the 
difference in track beginning and endpoint uncertainty. 

These selection criteriawere chosen to optimize the selection of cosmic ray muons with good point- 
ing resolution on the sky. The cut values were determined empirically using standard MINOS 
Monte Carlo events, inserting the moon, and maximizing the moon shadow. The number of muons 
that survived each cut is shown in Table[TJ Of the initial 83.54 million triggers, 62.5% survived all 
cuts, leaving 52.19 million muons. 

Table 1 

Fraction of events that survive each pointing cut 



Cut 


No. Remaining 


Fraction Remaining 


Total Tracks 


83.54xl0 6 


1.0 


1. Data Quality Cuts flU 


68.91 xlO 6 


0.825 


2. Number of Planes < 10 


62.06xl0 6 


0.743 


3. Track Length < 1.55 m 


60.10xl0 6 


0.730 


4. \~$vtx - ~Jend\ >0.021 


52.19xl0 6 


0.625 



There are 17,389 muons in a 2° half-angle cone centered on the moon and 16,41 1 muons in a 2° 
half-angle cone centered on the sun. The reason for this difference is that the detector is at a high 
latitude, thus the number of muons collected in winter near the sun's location will be fewer than 
the number collected near the moon. The data set includes five full yearly cycles plus five months. 
Those extra months came after the summer solstice, so the amount of time the sun spends above the 
horizon (and the angle of the sun above the horizon) continues to decrease as the period progresses. 



2.2 Geomagnetic Field 



The earth's magnetic field produces a relative east-west shift in the apparent arrival direction of 

cosmic ray primaries from the direction of the moon. Older experiments ( 11311411511 ) used the 

"impulse approximation" to calculate the expected lOTeV proton deflection, 0.15°. Since the path 

a particle travels is fairly long, the "impulse approximation" may not be valid over such large 

distances. The geomagnetic deflection Aa of a given cosmic ray with charge Z and momentum 

_± 

p(TeV/c), depends strongly on the local magnetic field conditions, B and the particle path / . It 
can be written: 

Z r d ™ —¥ —¥ 
tan(Aa) = — B x d I , (1) 
p Jo 

where d m = 384 x 10 3 km is the distance from the earth to the moon. To measure the effect of 
geomagnetic deflection on the moon shadow the magnetic field encountered by a cosmic ray pri- 
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mary between the moon and the earth was integrated along the cosmic ray path. The integral was 
calculated numerically, with the cross product computed at 20km intervals. Muons with arrival 
directions within a 2° half-angle cone of the moon's location were chosen from the data set. The 
geomagnetic field from the surface of the earth to 600 km was found using the International Ge- 
omagnetic Reference Field (IGRF) calculation [25]. The geomagnetic field from 600 km to the 
magnetopause (which varies depending on the IMF, here taken to be 70,000 km) was found using 
GEOPACK 11261271 . The geomagnetic field from the magnetopause to the location of the moon 
(384,000 km) was held constant at the value at 70,000km. The muons that survive to the depth of 
the far detector are usually the decay products of mesons created in the first cosmic ray interac- 
tion, typically at a height of 20 km. The integrals were thus calculated using ji + with momentum 
p^ = O.lpp from 0-20 km and protons from 20-384,000 km. The result of this calculation can be 
seen in Fig. [U where the surface energy of each muon was set to 0.7 TeV, the vertical muon thresh- 
old energy of the far detector, and that p p = lOp^. The maximum deflection was ARA = 0.02°, 




ARAcos(Dec )(°) 

Fig. 1. The distribution of integrated geomagnetic deflections for cosmic rays that pass near the moon in 
horizon coordinates. This distribution is for the case where the muon energy was determined from the far 
detector overburden profile and each primary cosmic proton had energy 10p M . 



and the mean deflection was ARA = 0.005°. Note that while there was a clear positive shift in RA, 
the distribution of deflections in Dec was more uniform. Any displacement of the moon shadow 
beyond 0.015° cannot be caused by the geomagnetic field, and it is unlikely that the geomagnetic 
field could cause an average deflection of greater than 0.01°. 
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2.3 Multiple Muons 



In addition to the smearing effects of the geomagnetic and interplanetary magnetic fields, the cos- 
mic muon angular resolution of an underground detector is limited by multiple Coulomb scattering 
in the surrounding rock as well as the geometric resolution of the apparatus itself. The geometric 
resolution of the detector is of order ±0.15° but the dominant effect 2070 mwe underground is mul- 
tiple Coulomb scattering. The angular resolution of the detector was measured using 3.12 million 
multiple muon events with reconstructed multiplicity equal to two, collected from August 1, 2003 
to December 31, 2007. Multiple muon events are created by a cosmic ray of sufficient energy to 
generate more than one energetic muon. The energy spectrum of multiple muon is similar to that 
of single muons, which is a steeply falling power law with a spectral index approximately mi- 
nus three. Most muons (multiple muon events included) are close to the low energy cutoff. The 
transverse momentum of these pairs is negligible compared to their longitudinal momentum, and 
they have angular separation <0.05° lfT4l at creation. The measured angular separation between 
these pairs quantifies the resolution of the detector. Pairs were selected by the same criteria used 
for single muons (see Sec. 12.11) . A total of 1.77 million multiple muon events survived. The dis- 
tribution of their angular separation ip, divided by \[2 to apply to the single muon resolution, is 
shown in Fig. [2l The resolution is defined to be the single muon angular separation within which 



a) 
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Fig. 2. a) The square of the single muon angular separation of multiple muon pairs, psi, in deg 2 . The shaded 
region represents 68% of the distribution, which is taken to be the resolution of the detector, ij)/\2 = 0.62°. 
The peak region of the distribution was fitted with a Gaussian function, shown by the dashed line. The fit 
parameter a = 0.322 ± 0.002°. b) The square of the single muon angular separation of multiple muon pairs 
out to 0.375 deg 2 . 

68% of the distribution lies. This value, shown by the shaded region in Fig. [2l is 0.62°. The peak 
region of the distribution was fitted with a Gaussian function, shown by the dashed line. The fit 
parameter a = 0.322 ± 0.002. The Gaussian function fits well over the first few bins, which de- 
scribes the bulk of the distribution, but the long tail begins to deviate from the Gaussian function at 
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ip 2 /2 ~ 0.06 deg 2 . The long tail can be attributed to Moliere scattering in the rock overburden lfi~3l . 



2.4 Monte Carlo Simulation 



The backgrounds for the shadow analyses were calculated using a Monte Carlo simulation. A 
muon arrival direction was chosen out of the known distribution of events in the detector (in hori- 
zon coordinates) and paired with a random time chosen from the known time distribution to find 
the muon's location in celestial coordinates. This was done for every data muon to create one 
background sample, and 1,000 background samples were created. 

A template of the expected distribution of the sun and moon shadows in celestial coordinates (right 
ascension (RA) and declination (Dec)) was simulated using the detector resolution determined 
from the multiple muon events in Sec. l2.3l to account for the smearing effects of coulomb scatter- 
ing and detector resolution. Cosmic rays were generated traveling towards a disk the size of the 
moon or sun. If the cosmic ray intercepted the disk it was removed. If not, the cosmic ray was 
assumed to produce a muon that would travel underground to the detector, and an angular devi- 
ation was selected at random from the multiple muon distribution. The resulting distribution in 
ARA ■ cos(Dec) (cos(Dec) normalizes for solid angle) and ADec can be seen for the moon in 

Rg. m 




Fig. 3. The expected cosmic ray shadow of the moon (from simulation) as seen in the far detector, in two 
different views: a) in three dimensions, and b) in a two dimensional projection of a). It is shown in celestial 
coordinates centered on the location of the moon. The color axis shows the surviving fraction of cosmic-ray 
induced muons. 
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3 Moon Shadow 



3.1 Two Dimensional Shadow 



A two dimensional maximum likelihood grid search was used to find the location of the moon. 
The muons in the background samples generated in Sec. 12.41 were sorted by separation from the 
moon's location in equal solid angle bins 0.01° on a side. The x coordinate is given by ARA = 
(RA^ — RAm) ■ cos(DecM), where the subscripts [i and M denote muon and moon, respectively, 
and the y coordinate is given by ADec = (Dec^ — Dccm)- ARA is modified by cos(D6Cm) to 
account for the projection of sphere onto a flat grid. One thousand cosmic ray background samples 
were averaged to create a smooth, isotropic background grid. The data were sorted in a similar grid. 
A grid search utilizing a log-likelihood method was employed to find the most probable position of 
a moon-like deficit. This method was invented by Cash ESI , first applied by COS-B 11291 , and first 



applied to the moon shadow by MACRO Ifl4ll . The moon shadow template (Sec. 12.41) was placed 
at a fixed position (x s ,y s ) on the data grid. The search was then performed by integrating over 
the the shadow template at this location. The template was then moved to adjacent bin location 
(x s ,y s +i), and the search was repeated. This process was repeated until the entire 4° by 4° data 
grid was scanned. The shadow that fit the data best was found by maximizing the shadow strength 
I s using the likelihood function: 



l(x,y,I.) = 2 y £ 

i=l 



ATobs 

nrth _ probs , jyobs , _J_ 



(2) 



where N- h = N^ ack — I s ■ P s (xi, yi) is the number of events expected in bin i, N^ ack is the number 
of muons from the smoothed background grid in bin i and P s (xi, yi) is the fraction of the cosmic 
rays at location (xi, yi) blocked by the moon, and is equal to one minus the distribution in Fig [3] 
I s ■ P s (xi, y^ is the number of events removed from bin i by the moon. To determine the strength 
of this deficit, the parameter A was defined as: 

A = l(x,y,0) - l(x,y,I 8 ), (3) 



which is a measure of the deviation from the null (no-moon) hypothesis. 

The two dimensional distribution of these deviations was drawn on a 4° x 4° grid, binned in 0.01° 
on a side, and can be seen in Fig. |4] in celestial coordinates. The greatest deficit is A = 30.9, 
centered on (-0.11 ± 0.09°, -0.13 ± 0.08°). This is consistent with the Soudan 2 fll3] shadow, 
which was offset by 0.1°E-W and 0.1°N-S. The error and maximal value was found by drawing a 
1 a contour around the bin location of A max . The value of I s at this location was 0.13, in which 
235.2 events were removed by the moon, which is consistent with the expectation of 297.8 events 
removed. In the high statistics limit, the distribution of A is the same as for a xl distribution ||2~8~1 . 
where v is the number of free parameters. In this case there is only one free parameter, I s , which 
means A = 30.9 has a significance of 5.6 a. 

A second method was used as a cross-check to find the significance of the observed deficit. Moon- 
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Moon 




ARAcosfDec )(°) 

M 

Fig. 4. The two dimensional moon induced muon deficit in 0.01 ° wide bins, in celestial coordinates. The x 
marks the expected location of the moon, and the dashed circle is the apparent size of the moon as viewed 
from earth. The cross marks the observed location of the moon, and the solid ellipse denotes the uncertainty. 
The greatest deficit is A = 30.9, a 5.6 a significance, centered on (—0.11°, —0.13°) 

less grids were created using the same background method described in Sec. 12.41 Searching through 
many moonless grids gave the probability of randomly finding a moon-like deficit for a particular 
value of A. Each 4° x 4° grid allowed 160,000 searches, for a total of 161 .4 million searches. There 
were no searches that had A > 23.0 and the distribution of A followed the expected \\ distribution. 



3.2 One Dimensional Shadow 



As a check on the resolution of the detector, a one dimensional moon shadow search was per- 
formed. The differential density of cosmic rays obscured by the moon, AN^/ACl, can be written 
as a two dimensional Gaussian convolution in polar coordinates (r,0): 

= X ( 1 - JL [ Rm r > dr ' /* d( f> e -(r> + r*-2rr> COsm A ^ 



AQ V 27TC7 2 Jo 
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where A is the average differential muon flux, a accounts for smearing from detector resolution, 
multiple Coulomb scattering and geomagnetic deflection, and R m =0.26°, the radius of the moon 
or sun. Performing the integration over gives: 



An 



A 1 - 



cr 2 Jo 



Rn 



r' dr'lo 



e 2" 



(5) 



where I is a modified Bessel function of the first kind. The integral that remains has no closed- 
form solution; however I has a rapidly decreasing Taylor expansion: 



lo(x) 



E 



>2/v 



2 4 fi 

H 1 1 

4 64 2304 



fc=0 4 fc (^!) 2 

Substituting the expansion and evaluating the integral gives 



+ 



(6) 



AN» 

An 



A 



1 _ Sn P -e 2 /2* 2 ( i _, 
2a 2 + 



? 2 - 2( x 2 X 



+ 



? 4 - 8#V 2 + 8a 4 ) m 



M 



192a 8 



(7) 



with = r since the observable in this case is an angular separation. This formula automatically 
produces a deficit of vri?^A events due to shadowing. The significance of the deficit can be found 
by fitting to Eq.[7Jand finding the difference between this x 2 value and the x 2 value obtained by a 
linear fit to the same data. 



3.3 Search in One Dimension 



The reconstructed muon angular separation from the moon or sun, AO, was binned in S^- m =0.10° 
increments. Since radial distance from the center is measured over a two dimensional projection, 
the solid angle of bin (i) increases when moving out from the center as AVti = (2i — 1) ■ S^ in 7r. 
Weighting the number of events in each bin by the reciprocal of the area resulted in the distribution 
Nil ' AVti, the differential muon density. 

Since the location of the moon was found to be offset from (0°, 0°), the arrival direction of each 
muon was adjusted by (0.11°, 0.13°) before performing the one dimensional moon shadow search. 
The position of the moon on the sky and the separation of each muon from the moon {AO), was 
found using the method described in Sec. 12.11 The AO distribution is shown in Fig. [5] with statistical 
error bars. There is a significant deficit in the location of the moon as AO — > attributed to the 
moon's blocking of the primary cosmic rays. A linear (no moon deficit) fit gives a y^/ndf = 
57.8/39. A fit to Eq.|7]( moon induced deficit) gives a x%l n< ^f — 25.3/38. The improvement in x 2 
of 32.6 corresponds to a significance of 5.6 a. The resolution found by this method is cr = 0.34 ± 
0.04°, which is consistent with a Gaussian fit to the multiple muon data shown in Fig. [2j When the 
same analysis is performed without shifting the muons by (0.11°, 0.13°), the significance of the 
moon shadow is reduced by 2.7 cr, and the shadow is more smeared out, with cr = 0.39 ± 0.06°. 
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Fig. 5. The differential muon flux with respect to displacement from the moon's location, binned in 
A9 = 0.1°. The dashed curve is the result of a linear (no moon effect) fit, which gives x^/ndf = 57.8/39. 
The solid curve is the best fit from Eq.|7] The Gaussian (moon-induced deficit) fit gives x^/ndf = 25.3/38, 
with parameters A = 1171.1 ± 4.8 and a = 0.34 ± 0.04°. 

4 Sun Shadow 



As viewed from earth, the sun also obscures a 0.52° diameter disk, approximately the same size 
as the moon. The one dimensional shadowing procedure from Sec. |3]and the two dimensional log- 
likelihood analysis described in Sec. 13. II were performed for cosmic ray muons coming from the 
direction of the sun. Historically, this has been a more difficult |fi~5l measurement to make because 
of the variability of the sun's magnetic field described in Sec. [5] 

The two dimensional sun shadow can be seen in Fig. [6] The sun A max = 14.6 occurs at (—0.29 ± 
0.13°, 0.27 ± 0.14°) and has a 3.8a significance. The value of I s at this location was 0.08, in 
which 144.6 events were removed by the sun, which is about half of the expected of 279.8 events 
removed. The sun shadow has a K max value that is less than half the A max value of the moon, and 
the shadow appears elongated along a line through the origin, rotated about halfway between the 
north and to the west. These features could be attributed to a smearing effect that is unobserved 
in the moon shadow, originating from the longer exposure to the IMF along the much greater path 
between the sun and the earth. 

The one dimensional sun shadow can be seen in Fig. |7J Note that the muons were shifted by 
(0.29°, —0.27°) to account for the offset from the expected location of the sun given by the two di- 
mensional sun shadow. There is a significant deficit in the location of the sun attributed to the sun's 
blocking of the primary cosmic rays. The improvement in x 2 of 15.1 (xi/ndf = 67.9/39, Xxs/ndf = 
52.8/38) corresponds to a 3.9 a significance . The Gaussian fit parameter a = 0.48 ± 0.07 for the 
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Fig. 6. The two dimensional sun induced muon deficit in 0.01 deg 2 bins, in ecliptic coordinates. The x 
marks the expected location of the sun, and the dashed circle is the apparent size of the sun as viewed from 
earth. The cross marks the observed location of the sun, and the solid ellipse denotes the uncertainty. The 
greatest deficit is A max = 14.6, a 3.8 a significance, centered on (—0.29°, 0.27°). The maximum of the 
color scale is set to A = 31 to allow for easy comparison with the moon shadow. 

sun shadow is somewhat larger than the value found by the fit to the moon shadow, a = 0.34±0.04. 



5 Interplanetary Magnetic Field Effects 

5. 1 Moon Shadow 

The IMF could have some effect on the moon shadow by the effects described in Sec. \T\ This 
could be observed by dividing the data into separate day/night samples as cosmic rays experience 
the maximal IMF difference when approaching the earth from the day and night sides. The effect 
of the IMF on the geomagnetic field is to produce a reduced geomagnetic field at night (12612713011 . 
This effect was observed by MACRO [fl5l . Daytime is defined here as when the sun's zenith angle 



14 



C\J 

'o>1100|r 
gpOOOb 

"D 

Z 900 N 



800:- 



700:- 



600]- 








t.MH+<' + tj 



_i i i i_ 



+ 



_i i i i i i i_ 



_i i i i_ 



3 4 

A6(deg) 



Fig. 7. The differential muon flux with respect to displacement from the sun's location, binned in 0.10°. A 
linear (no sun effect) fit gives xl/ nd f = 67.9/39. The solid curve is the best fit from Eq. |7] The Gaussian 
(sun-induced deficit) fit gives xh/ nd f = 52.8/38, with parameters A = 1040. 9±4.6 and a = 0.48±0.073°. 

was less than zero. The day-time sample contained 8,270 muons within a 2° half-angle cone around 
the moon, while the night-time sample contained 9,213 events. The reason for this difference is that 
the detector is only down for maintenance during the day, coupled with the fact that more of the 
data were collected in the fall and winter than spring and summer (see Sec. 12.11) . Integrating over 
five and a half years accounts for the 10% increase in muons collected near the moon at night. The 
moon shadow observed at night can be seen in Fig. [8] (Left), while the moon shadow during the day 
can be seen in Fig. [8] (Right). The center of the deficit is at (-0.09 ± 0.11°, -0.22 ± 0.13°) with 
^raax = 23.7 (4.9 a) for the data taken at night. For the data taken during the day, the center of the 
deficit is at (-0.21 ±0.12°, -0.07±0.12°), withA^ = 11.2 (3.3 a). A*g£* is consistent with the 
expected location of the moon in ARA, while the shadow during the day is shifted by 1.7 a. The 
shift in ARA at night is consistent with zero. The shift in ADec at night, however, is further from 
the expected location of the moon by 1.7 a. It is unlikely that the geomagnetic field could cause a 
deflection of greater than —0.02° ( Sec. 12.21) and the geomagnetic field should not cause a diurnal 
effect. A possible explanation is that the IMF has a general defocussing effect on cosmic rays, and 
cosmic rays detected at night experience less of this IMF effect because they arrive at earth from 
the opposite direction to the sun. 



5.2 Sun Shadow 



A greater IMF means a less prominent sun shadow, so the significance of the observable shadowing 
caused by the sun should decrease as the IMF increases. Since the IMF is caused by solar activity, 
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Fig. 8. The distribution of A values in celestial coordinates when the moon was visible during the 
night (Left), and when the moon was visible during the day (Right). The x marks the expected location 
of the moon, and the dashed circle is the apparent size of the moon as viewed from earth. The cross marks 
the observed location of the moon, and the solid ellipse denotes the uncertainty. The deficit has K max =23.1 
(4.9 cr) at night, and the center is at (-0.09 ± 0.11°, -0.22 ± 0.13°). The deficit has k max =\\2 (3.3 a) in 
the day, and the center is at (-0.21 ± 0.12°, -0.07 ± 0.12°). 

the significance of observable shadowing should increase as the number of sunspots, one measure 
of solar activity, decreases. The last maximum of the 1 1 year solar cycle occurred in 2001, and the 
next minimum occurred in December, 2008 [|23l . To search for a correlation between solar activity 
and strength of sun shadowing effects, the data were divided into five separate one dimensional 
grids of roughly equal statistics, and fit with both a linear function and a one-parameter Gaussian 
(Eq. [7]), holding the number of removed muons constant. The one dimensional search is used be- 
cause it requires fewer events to produce a statistically significant shadow and the sun was shown to 
be close to the center of the grid (see Fig. [6]). There were not enough accumulated muons to perform 
five of the two dimensional likelihood searches. The periods were Aug. 1, 2003 - Sep. 30, 2004, 
Oct. 1, 2004 - Oct. 31, 2005, Nov, 1, 2005 - Nov. 30, 2006, Dec. 1, 2006 - Dec. 31, 2007 and 
Jan. 1, 2008 - Dec. 31, 2008. These graphs can be seen in Fig.|9] 

The results of these fits are summarized in Table |2] There is a decrease in the size (a) and cor- 
responding error of the sun shadow as the sun approaches solar minimum. This decline could 
be correlated to the decrease in solar activity. Surprisingly, the period nearest to solar minimum, 
January 1, 2008-December 31, 2008, shows an increased size. For comparison, the yearly moon 
shadows for the same periods were calculated and summarized in Table |2] The moon shadow has 
more constant size over the five periods than does the sun. 

A more useful correlation between the sun shadow and the IMF is to perform two dimensional 
searches on data sets subdivided according to the vector components of the IMF, but there are 
problems with this approach. First of all, the number of data subsets is limited to two, because the 
2D log-likelihood analysis is statistically limited. Second, the IMF is known to have a complicated 
and time-varying structure. The ACE, Wind and IMP 8 spacecraft, which monitor the IMF while 



16 



8/01/03-8/31/04 




10/01/05-10/31/06 

350 



2 3 4 

A9(deg) 




2 3 4 

A0(deg) 

12/01/07-12/31/08 
ST* 350 n 



-g 300 
^ 250 



5- 200 
150 
100 
50, 



9/01/04-9/31/05 

«sr 350 P 
'o 

-g 300 



t3 250 



200 
150 
100 
50, 





11/01/06-11/30/07 

- 350 



2 3 4 

A9(deg) 




2 3 4 

A9(deg) 



'0 



A9(deg) 



Fig. 9. The differential muon flux with respect to displacement from the sun's location as a function of year, 
binned in 0.1°. The dashed curve is a linear (no sun effect) fit, while the solid curve is the best fit from Eq.[7j 



orbiting about the LI Lagrange point (225 earth radii in front of earth), report data averaged over 
one hour PDl . These data report the IMF in Geocentric Solar Ecliptic (GSE) coordinates, where 
the earth-sun line defines the x-axis, and the ecliptic north pole is the z-axis. It is difficult to 
adequately account for the variations in the field with these data, as was done in Sec. 12.21 for the 
geomagnetic field. The greatest difficulty is the fact that the IMF data is only available every hour 
for one location in space, whilst it is well known that the IMF changes somewhat rapidly in time 
and has a sectorized structure ETTl . The direction of the IMF is either "toward" or "away" from the 
earth, and the two field directions are separated by a thin layer known as the neutral current sheet, 
which passes the LI point very briefly as the x component of IMF reverses polarity. 
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Table 2 



The size of shadowing observed in each year's sun distribution. Also included is the corresponding moon 
shadow size. The significance is given in units of standard deviations (s.d.). 





Sun 


Moon 


Distribution 


a 


Xg 


Significance 


a 


Xg 


Significance 


Aug. 1,2003 -Sep. 30, 2004 


0.61 ± 0.28 


34 


1.6 s.d. 


0.34 ±0.11 


21 A 


1.9 s.d. 


Oct. 1,2004 -Oct. 31, 2005 


0.31 ±0.12 


42.3 


1.6 s.d. 


0.33 ± 0.08 


24 


4.0 s.d. 


Nov. 1, 2005 - Nov. 30, 2006 


0.27 ± 0.07 


45.5 


3.2 s.d. 


0.19 ±0.04 


38.6 


3.9 s.d. 


Dec. 1,2006 -Dec. 31,2007 


0.34 ±0.12 


42.3 


2.5 s.d. 


0.33 ±0.11 


32.2 


3.6 s.d. 


Jan. 1,2008 -Dec. 31,2008 


0.70 ± 0.29 


37.2 


2.1 s.d. 


0.26 ± 0.06 


42.6 


3.4 s.d. 



These IMF data were used to separate the muon data and search for IMF induced effects. The 
data were first separated into two sets by the magnitude in the expected deflection plane, which is 
perpendicular to the earth-sun line. The "large magnitude" set was where \B±\ > 4.2 nT; less than 
4.2 nT was the "small magnitude" data set. No significant translation of the shadow nor focusing 
effects were observed. The shadows had magnitudes equal to roughly half of the total sun shadow, 
consistent with what was expected from dividing the data set in half arbitrarily. 

An interesting effect was observed when the data were divided into sets where the IMF component 
along the earth-sun line (B x ) was greater and less than zero. For the B x > case, the most 
significant deficit was centered on (—0.07 ± 0.14°, 0.10 ± 0.12°), which is consistent with the 
expected location of the sun. This shadow shows none of the elongation seen in the total sun 
shadow (Fig. [6]). For the B x < case, the most significant deficit was centered on (—0.41 ± 
0.17°, 0.47 ± 0.20°). In fact, this shadow has almost no significance near the expected location of 
the sun, and has two deficits of almost equal significance an equal distance from the center. This 
offset is 2.4 a in RA and 2.3 a in Dec from the location of the sun. The sectorized structure of the 
IMF causes a cosmic ray traveling along the path from the sun to the earth to experience magnetic 
deflections in a multitude of directions. The rotating nature of the IMF adds an additional element 
to the magnetic deflection. It may be that the "inward" and "outward" separation of the cosmic ray 
data set along the earth-sun line is a more consistent predictor of cosmic ray deflection than the 
sparse measurements of the IMF in the y and z directions. 



6 Conclusions 

Using 52.3 million muons accumulated over 1857.91 live-days, the MINOS far detector has ob- 
served the cosmic ray shadow of both the moon and sun with high significance. The two dimen- 
sional moon shadow was seen with a significance of 5.6 a, centered on (—0.11 ± 0.09°, —0.13 ± 
0.08°), which suggests that the absolute pointing of the detector on the sky is known to 0.17±0.12°. 
The significance of the one dimensional moon shadow increased after shifting the muon arrival di- 
rections by (—0.11°, —0.13°) to be consistent with the observed location of the moon. The cosmic 



18 



Sun, IMF >0 Sun, IMF <0 

Y Y 




Fig. 10. The sun shadow for B x > (Left) and B x < (Right). The x marks the expected location of the 
sun, and the dashed circle is the apparent size of the sun as viewed from earth. The cross marks the observed 
location of the sun, and the solid ellipse denotes the uncertainty. The most significant deficit was located at 
(-0.07 ± 0.14°, 0.10 ± 0.12°) and had A max = 10.2 (3.2 a) for B x > 0. The most significant deficit was 
located at (-0.41 ± 0.17°, 0.47 ± 0.20°) and had A max = 9.2 (3.0 a) for B x < 0. 



ray shadow of the sun over the same time period was measured in two dimensions with a signifi- 
cance 3.8 a, centered on (-0.29 ± 0.13°, 0.27 ± 0.14°). 

Searches were made for possible dependencies of the moon and sun shadow on the IMF but no 
three sigma effects were observed. At lesser statistical significance, effects were observed in the 
difference of the moon shadow during the day and night and in the variation of the sun shadow 
with time. The shadow of the sun showed some indication (though less than three sigma) of vari- 
ation with the x component of the IMF, which could be interpreted as focusing (for B x > 0) and 
defocussing (for B x < 0). These latter observations of short term correlations with the IMF were 
made at the limit of the experiment's statistical power, and it will be interesting for other detectors 
to confirm these signals. 
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